Production Efficiency of Ultracold Feshbach Molecules in Bosonic and Fermionic 

Systems 
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We investigate the production efficiency of ultracold 
molecules in bosonic 85 Rb and fermionic 40 K when the mag- 
netic field is swept across a Feshbach resonance. For adiabatic 
sweeps of the magnetic field, the conversion efficiency of each 
species is solely determined by the phase space density of the 
atomic cloud, in contrast to a number of theoretical predic- 
tions. Our novel model for the adiabatic pairing process, de- 
veloped from general physical principles, accurately predicts 
the conversion efficiency for both ultracold gases of bosons and 
of fermions. In the non-adiabatic regime our measurements 
of the 85 Rb molecule conversion efficiency follow a Landau 
Zener model, with a conversion efficiency that is characterized 
by the density divided by the time derivative of the magnetic 
field. 

PACS numbers: 05.30.Jp, 03.75.Ss, 05.30.-d, 36.90.+f 



The production of ultracold diatomic molecules is an 
exciting area of research, with applications ranging from 
the search for the permanent electric dipole moment [1] 
to providing unique experimental access to the predicted 
BCS-BEC cross-over physics [2]. A widely used produc- 
tion technique involves the association of ultracold atoms 
into very weakly bound (~10 kHz binding energy) di- 
atomic molecules by applying a time varying magnetic 
field in the vicinity of a Feshbach resonance. By using a 
slow adiabatic sweep of the field through the resonance 
into the region where bound molecules exist [3-5], sam- 
ples of over 10 5 molecules at temperatures of a few tens 
of nK have been produced from both quantum degener- 
ate two-component Fermi gases [6-8] and atomic Bose- 
Einstein condensates [9-11]. Over the past two years, 
several experiments have probed the unique and exotic 
properties of these molecules. In [12] the first coherent 
superposition of two chemically different species, atoms 
and Feshbach molecules, was demonstrated, and the van- 
ishingly small molecular binding energy was measured. 
Bose-Einstein condensation (BEC) has been observed in 
Feshbach molecules formed from two fermionic atoms 
[13,14], and these molecules provide a unique tool for the 
controlled investigation of the basic pairing phenomena 
relevant to superconductivity [15]. 

Despite the widespread use of this technique, the 
molecule production process has received little experi- 
mental attention, and a microscopic description that ac- 
curately predicts the conversion efficiency as a function 
of sweep rate, atom type, density and temperature has 



not yet been developed. Theoretical work on bosonic 
systems always assumes the existence of a condensate 
and hence 100% molecular conversion for sufficiently slow 
sweeps, although lower values are always observed in ex- 
periments due to finite molecular lifetimes. Meanwhile in 
quantum degenerate fermionic systems, conversion lim- 
its of 50% have been suggested in certain limiting cases 
[16] but these theories do not cover the full range of ex- 
perimentally accessible parameters [13]. In this paper 
we present conversion data in the adiabatic regime from 
both an ultracold but non-condensed cloud of bosonic 
85 Rb atoms and an ultracold cloud of 40 K fermions, at a 
range of degeneracies. We show that a single, very gen- 
eral description of the pairing process accurately predicts 
the relationship between conversion efficiency and phase 
space density in both cases. This theory demonstrates 
that the complete conversion of the condensate is just 
the limiting behaviour of a bosonic thermal cloud as the 
phase space density increases. 

Outside the adiabatic regime, several predictions have 
been made for conversion efficiency as a function of mag- 
netic field sweep rate [3] (and refs. therein) but have 
only been tested against experimental data over a nar- 
row range of parameters [6,7]. In the first section of this 
paper, we present a thorough investigation of the non- 
adiabatic regime using the 85 Rb system. Our results sup- 
port a general Landau Zener type theory for conversion 
as a function of magnetic sweep rate and density that is 
discussed in detail in [3,4]. 

Detailed descriptions of the 85 Rb experimental ap- 
paratus and the magnetic field sweeps for producing 
molecules are contained in [17] and [18], respectively. 
We use evaporative cooling to produce thermal clouds 
of 50,000-130,000 atoms at temperatures of 26-94 nK. 
The atoms are held in a purely magnetic 'baseball' trap 
[17], with fixed trapping frequencies of 17.6 x 17.6 x 
6.8 Hz. For efficient evaporation, the bias field is held 
at 162 G, where the scattering length a is positive. For 
slow magnetic field ramps, Rb2 molecules are only pro- 
duced when the field is ramped upward through the res- 
onance, which is located at 155 G; hence the first step 
in molecule production is to rapidly jump the magnetic 
field from 162 to 147.5 G. We then sweep the field back 
up to 162 G at a chosen linear rate, producing molecules 
as we pass through the Feshbach resonance. The cloud 
is then held at 162 G for 10 ms, during which time 
all the molecules undergo spontaneous one-body disso- 
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ciation and the two dissociated atoms rapidly leave the 
cloud with several mK of kinetic energy [18,19]. Thus the 
number of molecules formed is just the difference in atom 
number before and after the molecular conversion sweep. 
Given the unusually low densities of the cloud in our 
experiments (~10 n cm -3 ), two and three-body atomic 
decay and collisional molecular decay rates are negligi- 
ble and do not affect our measurement of the molecular 
fraction [18]. In other experiments, operating at higher 
densities (e.g. [10]), collisional decay cannot be ignored. 
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FIG. 1. (a) Molecular conversion in 5 Rb as a function 
of inverse magnetic field sweep rate. The initial conditions 
of the atomic cloud were N = 100, 000, T = 30 nK and 
n = 1.5 x 10 11 cm" 3 . The solid line is a fit to Eq.(l). (b) 
Initial mean density of the atomic cloud (n) divided by 1/e 
sweep rate (-Bi/ e ) and plotted against mean density. The 
dashed line shows the best fit of Eq.(2) to the data. For the 
four low density points, the vertical error bars are dominated 
by the uncertainty in fitting B 1 / e . At higher densities, the 
cloud experienced significant heating during the ramps across 
the resonance; hence the density of the final point has signifi- 
cant uncertainty. This heating limited the maximum density 
used. 

First we investigated the molecular conversion effi- 
ciency as a function of magnetic field ramp rate. A typ- 
ical data set is shown in Fig. 1(a) and is well fitted by a 
Landau Zener formula for the transition probability at a 
two-level crossing 



where N max is the asymptotic number of molecules cre- 
ated for a very slow ramp, B is the magnetic field sweep 
rate and /3 is a fitting parameter. j3/B is often called the 
Landau Zener parameter 5i z . 

The dependence of 5i z on the mean density of the cloud 
n, the magnetic field width of the resonance A, the back- 
ground scattering length a^g, and B can be derived intu- 
itively up to a constant a by considering the time taken 
to cross the resonance divided by the mean-field coupling 
strength, 



;Aa bg / B. 



(2) 



Equation 2 shows that for true Landau Zener be- 
haviour, the quantity n/Bx/ e should be constant, inde- 
pendent of density. (At B 1/e , S tz = 1 and N mo i/N rnax = 
63%.) Similar datasets to Fig. 1(a) were taken over a 
range of initial conditions and n/ Bi/ e is displayed versus 
n in Fig. 1(b). The data support a constant value for 
n/Bx/ e and hence a Landau Zener dependence on sweep 
rate. The value for a, extracted from the data of Fig. 1(b) 
is 4.5(4) x 10 -7 m 2 s -1 . In Ref. [3], the Landau Zener 
behaviour is rigourously derived for a zero temperature 
condensate. The authors predict a to be roughly 1/8 of 
the value that we extract from our thermal cloud data. 
This difference is currently being investigated [20]. 
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FIG. 2. Molecule conversion efficiency in 85 Rb as a function 
of peak phase space density (circles) . The solid line shows a 
simulation based on our conversion theory, fitted to the data 
with a single pairing parameter 71, = 0.44(3) (Eq.(3)). The 
dashed lines indicate the uncertainty in 75. Inset: Molecule 
conversion efficiency in 40 K as a function of peak phase space 
density (for full data range see Fig. 3) . The best-fit simulation 
based on our conversion theory is shown as a solid line, with 
the uncertainty indicated by dotted lines. 

Our second experiment used both bosonic and 
fermionic systems and investigated molecular conversion 
efficiency at a sweep rate that is slow (adiabatic) with 
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respect to B\/ e (e.g. slower than 150 /xs/G for the con- 
ditions of Fig. 1(a)) but still fast compared to the rate of 
molecule production via three-body inelastic collisions. 
This regime [21] is important for producing molecules 
with a limited lifetime. The conversion efficiency for 85 Rb 
is shown in Fig. 2 as a function of the peak phase space 
density of the cloud. Note that for adjacent points on 
Fig. 2 (i.e. points with very similar phase space den- 
sity), the number of atoms often varies by over a fac- 
tor of two, indicating that the conversion efficiency is 
related to phase space density but not to number or tem- 
perature individually. The inset plot in Fig. 2 shows 
that the same relationship between conversion efficiency 
and phase space density may be observed in ultracold 
fermionic 40 K, as discussed later in this paper. 

One can argue intuitively that there should be a close 
relationship between phase space density and conversion 
efficiency via the adiabatic sweep process. An adiabatic 
process smoothly alters the wavefunction describing atom 
pairs but does not change the occupation of states in 
phase space. Thus to form a molecule, a pair of atoms 
must initially be sufficiently close in phase space (small 
relative momentum compared to h divided by separation) 
that their combined wavefunction can evolve smoothly 
into the highest bound molecular state as the resonance 
is crossed. This theory predicts 100% conversion in a 
condensate because every atom occupies the same state 
in phase space and so is able to form a molecule with any 
other atom. For a non-condensed sample, we are still able 
to convert a substantial fraction to molecules because 
there are no constraints on the center of mass velocity of 
the atom pair; it simply converts to the velocity of the 
molecule. 

An analytical relationship between phase space den- 
sity and conversion efficiency is only possible in the limit 
where the conversion fraction is very small. For almost 
any finite efficiency, one cannot ignore that the forma- 
tion of molecules reduces the pairing options for those 
atoms that remain, leading to the saturation behaviour 
observed in the experimental data of Fig. 2. To account 
for this, we must add a monogamy clause to our model; 
two atoms will form a molecule if they satisfy the phase 
space proximity condition described above and neither 
has already paired with another atom to form a molecule. 
To test this model for the pairing process we used a 
Monte Carlo type simulation based on these two prin- 
ciples. For each run, we randomly assigned positions and 
momenta to typically 10,000 atoms [22] ensuring that 
the ensemble retained the correct phase space distribu- 
tion for the specified temperature and trap frequencies. 
The simulation program searched for a partner for each 
atom in turn, removing both atoms from future searches 
if the pairing condition was satisfied. The simulation 
ended when no more molecules could be formed. Atoms 
were considered sufficiently close in phase space to form 
a molecule if they satisfied the relation 



\Sr re i mSv re i\ < 7/1, (3) 

where Sr re i is the separation of the pair, m is the atomic 
mass, Sv re i is the relative velocity between the two atoms, 
and 7 is a constant to be determined by fitting the out- 
put of the simulation to experimental data. Finally, for 
the boson case, the simulated conversion fraction was in- 
creased by 2% to allow for two-body correlations that 
assist molecule formation in a cloud of identical bosons 
[23]. 

For each value of 7, the simulation was run over a 
range of phase space densities and compared to the ex- 
perimental data using a least squares routine. The best 
fit for the bosonic data was given by jb — 0.44(3) and 
is shown as the solid curve on Fig. 2. The error on this 
result is dominated by a 12% systematic uncertainty in 
the phase space density (resulting from uncertainties in 
atom number and temperature). The dotted lines indi- 
cate the range of values for the pairing parameter 75 that 
result from this phase space density uncertainty. The 
excellent agreement between experiment and simulation 
over a wide range of phase space density provides strong 
evidence in support of our molecular conversion model. 

Our model for the pairing process is based on very gen- 
eral arguments about adiabaticity and atom availability. 
With appropriate phase space distributions, it is equally 
applicable to molecule formation in a fermionic system. 
This provides another stringent test of its validity. To 
study the molecule conversion efficiency in a Fermi sys- 
tem, we create an ultracold, two-component Fermi gas 
as outlined in [24]. To achieve a wide range of phase 
space density we follow cooling of the gas in an optical 
dipole trap with a recompression of the trap and con- 
trolled parametric heating through a modulation of the 
trap strength. This allows creation of Fermi gases rang- 
ing in temperature from < 0.05 Tp to 1.3 Tp, where Tp 
is the Fermi temperature. The dipole trap that confines 
the final Fermi gas is characterized by radial frequencies, 
v r , between 312 Hz and 630 Hz and an aspect ratio of 
Vrlv z — 70±20. With this ultracold Fermi gas, we create 
molecules as described in [6] using a broad s-wave Fes- 
hbach resonance between the /, m/) = 1 9/2, —9/2) and 
|9/2, -7/2) states of 40 K located at B = 202.1 ± 0.1 G 
[15]. The magnetic field sweep used to create molecules 
starts at 202.83 G (where a < 0), ends between 201.10 G 
and 201.59 G and occurs at inverse sweep speeds ranging 
from 640 to 2900 /xs/G. Within this range the conversion 
fraction is independent of sweep rate, confirming that 
the conversion is adiabatic. On the time scale of these 
magnetic field sweeps, the loss of molecules and atoms is 
negligible [25]. At the end of the magnetic field sweep, 
we immediately turn off the dipole trap and allow the 
gas to ballistically expand. We then use one of two tech- 
niques to probe the molecular conversion: (1) We com- 
pare the number of atoms remaining on the a > side 
of the resonance, to the number of atoms measured at 
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a < as described in [6]. (2) We convert the remaining 
atoms in the |9/2, -7/2) state to the |9/2, —5/2) state 
by applying a radio frequency it pulse. We then disso- 
ciate the molecules by ramping the magnetic field back 
to the a < side of the Feshbach resonance. Finally, we 
image the |9/2,— 7/2) and 9/2, — 5/2) states separately 
to extract the molecule and atom numbers, respectively. 
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FIG. 3. Molecule conversion efficiency in 40 K as a function 
of T/Tf- (Note that T/Tf uniquely determines the phase 
space density of a fermionic cloud). The initial clouds have 
mean densities for each spin state ranging from 1 x 10 12 to 
2 x 10 13 cm~ 3 and 52(2)% of the atoms are in the m f = -7/2 
state. The data are fitted with the same conversion model 
as for 85 Rb. The best fit curve has a pairing parameter of 
0.38(4). The dotted lines indicate the uncertainty on this 
result. Note that conversion efficiencies far greater than 50% 
are measured. 

Experimental data for molecule formation efficiency in 
the fermionic gas is shown in Fig. 3, and a limited range 
of the same data is shown in the inset of Fig. 2. The 
data was simulated with the same molecular conversion 
model as the bosonic data, modified to provide the cor- 
rect Fermi-Dirac phase space distributions for a 50/50 
spin mixture and allow pairing only between atoms with 
unlike spins. The pairing parameter for the best fit sim- 
ulation is 7/ = 0.38(4), where the error includes a 10% 
systematic uncertainity in T/Tp. This result is in ex- 
cellent agreement with the value of = 0.44(3) from 
the bosonic system, especially given that both results 
are dominated by systematic rather than statistical un- 
certainties. This agreement provides clear evidence that 
our model and the physics underlying it accurately de- 
scribes molecular conversion efficiency in both bosonic 
and fermionic systems. One fundamental difference be- 
tween our model and many previous theories is that we 
allow each atom to pair with any other available atom 
that satisfies the phase space proximity condition. Other 
theories have considered only one potential partner for 



each atom, which has led to a ceiling for the fermionic 
conversion of 50% [16], contradicting the experimental 
data of Fig. 3. 

In summary, we have performed the first thorough ex- 
perimental investigation of the efficiency of molecular for- 
mation when the magnetic field is swept across a Fesh- 
bach resonance. In the non-adiabatic regime, we have 
shown that the conversion efficiency follows a dependence 
on density and sweep rate that is well approximated by a 
two level Landau Zener model. In the adiabatic regime, 
using an ultracold but non-condensed gas of bosonic 85 Rb 
atoms and a quantum degenerate gas of fermionic 40 K, 
we find that the conversion efficiency is monotonically 
related to the initial phase space density of the atomic 
cloud. We have shown that the pairing process allows 
each atom to form a molecule with any other atom in the 
cloud provided that they are sufficiently close together in 
phase space. The agreement of experimental data from 
both bosonic and fermionic systems with a single simula- 
tion provides compelling evidence for the validity of this 
molecule conversion model. 
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